General Relation between Entanglement and Fluctuations in One Dimension 
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In one dimension very general results from conformal field theory and exact calculations for 
certain quantum spin systems have established universal scaling properties of the entanglement 
entropy between two parts of a critical system. Using both analytical and numerical methods, we 
show that if particle number or spin is conserved, fluctuations in a subsystem obey identical scaling 
as a function of subsystem size, suggesting that fluctuations are a useful quantity for determining 
the scaling of entanglement, especially in higher dimensions. We investigate the effects of boundaries 
and subleading corrections for critical spin and bosonic chains. 
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Entanglement entropy, which measures nonlocal cor- 
relations in a quantum system, plays an important role 
in such diverse areas as the study of black holes [1] and 
quantum computation [2] . More recently, much attention 
has been focused on entanglement in condensed matter 
systems [3] and in particular on the role of entanglement 
in quantum phase transitions at zero temperature [4] . A 
significant discovery arising from this investigation has 
been the universal scaling of entanglement entropy in 
one-dimensional quantum critical systems described by 
conformal field theory (CFT) [1, 5-7]. Despite these ad- 
vances, the experimental relevance of these theories has 
remained unclear; the same feature that makes this quan- 
tity so universal — mainly, the fact that entanglement en- 
tropy is defined without reference to the observables of a 
system — has precluded its measurement in real quantum 
many-body systems. Recently, however, it was shown [8] 
that for the special case of free fermions the entanglement 
entropy can be related exactly to the full set of cumulants 
of the charge fluctuations, suggesting that entanglement 
entropy could be accessed through the fluctuations. 

In this Letter we propose that the fluctuations of a 
conserved charge is an interesting quantity to study in 
relation to entanglement entropy beyond the free-fermion 
case. In particular, we show that in one-dimensional crit- 
ical systems with particle number or spin conservation, 
the variance of the fluctuations in a subsystem (hence- 
forth simply "fluctuations") behaves very similarly to 
the entanglement entropy even when there are interac- 
tions. Fluctuations, like entanglement entropy, diverge 
logarithmically as a function of subsystem size in confor- 
mally invariant systems with a globally conserved charge. 

The entanglement entropy of a subsystem A of size x 
embedded in a larger system of size L is given by the von 
Neumann entropy S(x, L) = — Tr pA of the reduced 
density matrix pa for subsystem A. For a critical system 
described by a CFT with central charge c, the entangle- 
ment entropy at zero temperature for L — > oo was shown 
to have the universal behavior [5] 

S(x) = - lnx + sx, (1) 



where s± is a non-universal constant. For later compari- 
son it is useful to note that this was achieved by comput- 
ing the quantity Tr jy\ and differentiating with respect 
to n at n = 1. The versatility of this formalism lay in 
the fact that Tr p n A transforms simply under conformal 
mappings, allowing S(x,L) to be computed for finite L, 
finite temperature, and different boundary conditions. 

Now for the same setup, consider the number fluctua- 
tions in subsystem A, 

Fa = ((Na-(N a )) 2 ), (2) 

where Na counts the number of particles in subsystem A. 
For spins we replace N by S z . We expect that the fluctu- 
ations will behave similarly to the entanglement entropy, 
since for a pure state J- a = J~b where B is the remain- 
der of the system. As noted for entanglement entropy, 
this symmetry implies that generically, the fluctuations 
reside mainly on the boundary between the two subsys- 
tems, leading to an area law with possible logarithmic 
corrections [9]. Indeed, it is easy to show that all even- 
order cumulants of the particle number satisfy this prop- 
erty, which explains the absence of odd cumulants in the 
formula of Ref. 8. Moreover, Ta — for separable states, 
and for valence bond states J- 'a coincides (up to a factor 
1/4) with both the von Neumann and valence bond en- 
tropies [10]. On the practical side, for most systems in 
any dimension the fluctuations are easier to compute nu- 
merically than the valence bond entanglement entropy, 
which was introduced partly for its computational con- 
venience relative to the von Neumann entropy. 

We first consider Luttinger liquids (LLs), which de- 
scribe the low-energy properties of many one-dimensional 
systems [11]. From LL theory in the limit L — ¥ oo we 
have 7t 2 J-ll = ([4>( x ) ~ 0(O)] 2 ) where <fi is the "charge" 
field, so that at zero temperature 

n 2 F hh (x)=Kln^, (3) 

with K the Luttinger parameter and a a short-distance 
cutoff. As for entanglement entropy [5], the same result 
with K — > K/2 and x — > 2x is obtained when there is 
a boundary, due to the constraint </>(0) ~ constant. An 
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interesting confirmation of the LL result comes from the 
v = 1/m fractional quantum Hall states [12], for which 
K = v. A detailed calculation [13] shows that in the 
time domain, charge fluctuations across a quantum point 
contact with quantum Hall wires are given by Ti 2 T(t) = 
v\n{t/8) with short-time cutoff 5, as might be expected 
from Eq. (3) and Lorentz-invariance. 

The logarithmic scaling of fluctuations extends beyond 
LLs, and holds generally for critical models described by 
a CFT with a conserved U(l)-current (i.e., fixed total 
particle number or spin component), which is always de- 
scribed by a massless free boson. As for entanglement 
entropy these results therefore extend simply to finite 
size, finite temperature, and different boundary condi- 
tions via conformal mapping. To see this, note that if we 
define the characteristic function 



Ma(A) = ( e 'M^-<^») 



-ffA 2 /(2vr 2 ) 



(4) 



which transforms simply under conformal mappings, 
then k 2 T a = -n 2 M'Jl(0) = g\a{x/a). The prefactor g 
can always be fixed by considering the physical meaning 
of the charge, but we can give a heuristic argument for 
its value as follows. At finite temperature 1/(3 (we set 
h = ks = v = 1, where v is the effective velocity) the 
mapping z — > z' = ((3/2tt) Inz in Eq. (4) gives 



w 2 T(x,(3) = 5 In 
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For sufficiently large x ^> (3 such that interactions across 
the boundary can be neglected (which is possible since 
correlations decay exponentially), we may consider the 
subsystem A to be a grand canonical ensemble in equi- 
librium with a bath consisting of the remainder of the 
system [14]. This is of course only possible if the total 
particle number is fixed. Then from standard statistical 
mechanics one has F{x,(3) <~ kx/(3 where n = dn/dji is 
the compressibility (susceptibility \ = dm/dB for spins), 
so that by matching Eq. (5) for x » j3, a we find 
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(6) 



We have put in the velocity v for completeness. Note that 
this is consistent with the LL expression since K = irvn 
[11]. Eq. (6) turns out to be quite general, as we will see 
later. These results generalize the logarithmic scaling 
of fluctuations noted for noninteracting fermions [15] to 
conformally invariant, interacting systems. Interestingly, 
to leading order the entanglement entropy and fluctua- 
tions both obey logarithmic scaling, and 



S(x) c 



x^> a. 



(7) 



This generalizes Ref. 8, with the ratio modified by the 
central charge and compressibility, for the free bosonic 



theory where only the second cumulant is non-zero. Of 
course, there are subleading corrections to both quanti- 
ties, which we study in detail. In many cases the sublead- 
ing terms of the entanglement entropy and fluctuations 
also behave similarly, especially when there are bound- 
aries. Indeed, in the presence of boundaries the loga- 
rithmic prefactors for both quantities are half of their 
periodic values, so that Eq. (7) remains unchanged. 

In the following we study the detailed behavior of fluc- 
tuations in several important models, including systems 
that are not described by LLs. 

XXZ Model.— We first consider the spin-1/2 XXZ 
Hamiltonian 

tfxxz - £(Sf + SfSf +1 + AS?S? +1 ) (8) 

i 

for — 1 < A < 1 where the model is gapless (the point 
A = — 1 is not conformally invariant). At A = it 
reduces to an exactly solvable problem of free fermions. 
Using the Jordan- Wigner transformation we can use the 
methods of Ref. 16 to compute the correlation matrix 
Gij , which in the limit L — > oo is found to be G? bc = 6ij — 
sinc[7r(i — j)/2] for periodic boundary conditions (PBCs) 
and G° bc = % — sinc[7r(i— j)/2]+sinc[7r(z+j)/2] for open 
boundary conditions (OBCs), where sine x = 1 for x = 
and sinx/x otherwise. For exact diagonalization the 
slightly more complicated finite-L expressions were used. 
The entanglement entropy was numerically computed as 
in Ref. 6, while JxxM = £fj=i[<^> - (S?)(S*)] = 

( £ - EL=i Gij)/ 4 f or a block of £ sites. 

For I ^> 1, an analytical result was obtained for the 
entanglement entropy for PBCs [17]: 



Sxx.(t,L) = -log 2 £ + si, 



(9) 



where c = 1, s\ ~ 1.047, and we use log 2 for the entropy. 
For easier comparison to numerical data we always work 
with the formula for finite L which corresponds to the 
mapping I — > (L/ir)sm(ir£/L). Similarly, we find that 
the spin fluctuations for PBCs are given by 



Tr 2 T xx (£ 1 L) = lii£ + f 1 



(10) 



plus 0(£r 2 ) corrections, where j\ = 1 + 7 + In 2 and 
7 is Euler's constant. This is consistent with K = 1 
in the corresponding LL description. Eq. (10) was also 
derived in a different context [18]. Fig. 1 compares the 
exact diagonalization result to the analytical results for 
both the entanglement entropy and spin fluctuations for 
PBCs; even for 100 sites the agreement is excellent. 
For OBCs the spin fluctuations are given by 



1 



Fxx(2£,L) 



+ 



1 (-1)' 
2tt 2 (2^) ir 2 {2£) 



[\n(2£) + 7 + ln2] 



(11) 
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FIG. 1: (color online) Results from exact diagonalization for 
the entanglement entropy (squares) and spin fluctuations (cir- 
cles) for a spin- 1/2 XX chain with PBCs (upper symbols) and 
OBCs (lower symbols, inset shows zoom around the chain cen- 
ter), L = 100. For PBCs solid lines are the analytical results 
in Eq. (9) and Eq. (10). For OBCs, solid lines are the ana- 
lytical result (11) for the fluctuations and a fit to Eq. (12) for 
the entropy with c ~ 0.997 and si ~ 1.050 (the exact values 
are c = 1 and si ~ 1.047). The fluctuations are scaled by 4. 

plus 0(£~ 2 ) corrections. The result is very similar to the 
oscillating form found in Ref. 19 for the entanglement 
entropy in the presence of a boundary, with an addi- 
tional oscillating contribution oc (— 1) (\n£)/£. As shown 
in Fig. 1, the entanglement entropy is described well by 

S&UL) = §1^(20 + S i +ai ^- a 2 t±t. (12) 

The XXZ model is solvable for all A by Bethe ansatz, 
but extracting the exact fluctuations is not practical. 
However, from LL theory we know that for |A| < 1 its 
asymptotic behavior (for PBCs) is [11] 

ioz qz\_(qz \/<?*\ = — 1 i 2A ~ 2 ( ~ 1 ) r 

v- f i+r k - J i I N^i+J"/ \ i / 2^ 2 2 2K ' ' ' ' 

(13) 

where only the leading terms are shown with K = 
(1/2) [1 — (cos -1 A)/7r] _1 and A 2 a non-universal con- 
stant. From Eq. (13) we find 

(-lY 

n 2 Txxz{e)=K\n£ + f 2 -A 2 ^^- (14) 

plus 0(£~ 2 ) corrections, while the entanglement entropy 
is given by Eq. (9) with c = 1 but a different constant. 
In the derivation of Eq. (14) a term proportional to £ was 
suppressed, since it arises from the short-distance physics 
not taken into account by Eq. (13) and we are guaranteed 
by Eq. (3) that the leading term is oc ln£. It is interest- 
ing that the logarithmic divergence arises from the 1/r 2 
term in the correlation function, which for K < 1 is the 
subleading contribution at large r. Thus the diverging 
fluctuations are due to short-distance correlations. 
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FIG. 2: (color online) DMRG results (circles) for the spin 
fluctuations in a spin-1/2 XXZ chain with PBCs, L = 100. 
Only A = —0.2 to 0.9 are shown (from top to bottom, in 0.1 
increments). Solid lines are fits to Eq. (14). (INSET) Fitted 
Luttinger parameter K for |A| < 0.9 (dots), with the solid 
line showing the Bethe-ansatz solution. The first and last 10 
sites were dropped for fitting purposes. 

Since K = 1 for the XX model the oscillating term 
can be neglected to 0(£~ 2 ) in agreement with our pre- 
vious result (10). The same is true for A < 0, because 
K > 1. In contrast, for A > the oscillations grow 
larger as we approach the Heisenberg point A = 1 where 
K = 1/2. Interestingly, although the same oscillating 
terms are not present in the entanglement entropy itself, 
they were recently shown to be a feature of the Renyi en- 
tropies [20]. At exactly the Heisenberg point the other- 
wise irrelevant umklapp term cos(4</>) becomes marginal, 
and the oscillating part of the spin-spin correlation func- 
tion acquires a logarithmic correction A' 2 {— 1)' 'y/hir/r. 
In this case the corresponding term in the fluctuations is 
oc (— l) £ \/ln £/£. More importantly for finite L the Lut- 
tinger parameter K is also renormalized. Fig. 2 shows 
the fit of density-matrix renormalization group (DMRG) 
[21] data to Eq. (14), with excellent agreement between 
the fitted K and the Bethe ansatz solution for |A| < 0.9. 

Other spin chains. — The isotropic Heisenberg model 
at A = 1 is also an example of the SU(2) Wess- 
Zumino-Witten (WZW) nonlinear cr-model with (inte- 
ger) topological coupling constant k [22], which de- 
scribes the low-energy physics of many other spin chains. 
Critical Heisenberg chains with half-integer spin be- 
long to the k = 1 universality class with central 
charge c = 1 [23], for example, and we can deduce 
tt 2 J-wzw{x) ~ (fc/2) ln(x/a). An interesting case where 
we can compute the fluctuations analytically is the 
Haldane-Shastry (HS) model with 1/r 2 interactions [24], 
which has exactly known spin-spin correlation function 
(§Ur%)-(%+r)(§t) = (-l) r Si(7rr)/(47rr) with Si(x) = 
J Q dt (sin t) 1 1. For large I we find 

(15) 
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plus 0(£~ 2 ) corrections, where is an integral whose 
value is / 3 ~ 0.197. This is consistent with the WZW 
fixed-point predictions, and in particular with the spin- 
1/2 Hcisenberg chain without the umklapp term. This 
also shows that the logarithmic scaling is not affected by 
long-range interactions. 

A case where k > 1 is the spin-s Takhtajan-Babujan 
(TB) chain [23]. Using DMRG we have checked the case 
k = 2s = 2 and hence 7t 2 J : "tb(^) ~ ln(x/a), with cor- 
responding central charge c = 3k/ (2 + k) = 3/2. This 
also explicitly confirms the result g — irvx — s = fc/2 
in Eq. (6). Another interesting example is the open- 
boundary Uimin-Suthcrland model which is a critical 
spin-1 chain with SU(3) symmetry because both the en- 
tanglement entropy and spin fluctuations exhibit oscilla- 
tions with a period of 3 sites due to the higher symmetry. 

Bose-Hubbard Model. — The Hamiltonian for the one- 
dimensional Bosc-Hubbard model is [25] 

#BH - -tJ2(blb i+1 + h.C.) + j ^hiihi ~ 1), (16) 



where bi is the bosonic annihilation operator on site i, 
hi = b\b i , t is the tunneling amplitude and U is the on- 
site repulsion. The model describes interacting bosons 
on a lattice, and can be realized with cold atoms trapped 
in an optical lattice [26]. As an example, we consider the 
case of half-filling since it approaches the XX model as 
U — > oo. The low-energy physics of the model is that of 
a LL [11], so that the density fluctuations are given also 
by Eq. (14). However, because K > 1 the oscillations 
are absent for PBCs and very weak for OBCs. We have 
confirmed all of these results with DMRG and obtained 
curves similar to those shown in Fig. 1 for the XX model. 

Gapped models. — For gapped models, we expect that 
the fluctuations will obey an area law like entanglement 
entropy, although the ratio is no longer fixed by con- 
formal arguments. This can be checked explicitly for 
the Affleck-Kennedy-Lieb-Tasaki model [27], for exam- 
ple, where analytical results are available but also the 
valence bond picture makes the relation intuitive [7]. 

Conclusion. — From the theory of quantum critical phe- 
nomena one expects an intimate relation between entan- 
glement and fluctuations. In one dimension this expec- 
tation is borne out by Eq. (7): like entanglement en- 
tropy, the number fluctuations scale logarithmically in 
critical models described by a CFT with a globally con- 
served charge. Moreover, our detailed investigation of 
the effects of boundaries and subleading corrections for 
several important models suggests that studying fluctua- 
tions is a powerful approach to understanding the scaling 
of entanglement entropy more generally. This has clear 
advantages: First, fluctuations are accessible in experi- 
ments, perhaps most easily for cold atoms in an optical 
lattice. Second, from the computational point of view 
they are easier to calculate, both analytically and nu- 



merically, than entanglement entropy. In particular, the 
close connection between entanglement entropy and fluc- 
tuations suggests that investigations of the latter with 
quantum Monte Carlo for d > 1 will provide clues to 
such important questions as whether the L d ~ Y InL scal- 
ing of both S and T in free fermions [15] holds more 
generally for interacting fermions [28]. 
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